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Abstract 

We derive a nonlinear integral equation (NLIE) for some bulk excited states of the 
sine-Gordon model on a finite interval with general integrable boundary interactions, 
including boundary terms proportional to the first time derivative of the field. We 
use this NLIE to compute numerically the dimensions of these states as a function of 
scale, and check the UV and IR limits analytically. We also find further support for 
the ground-state NLIE by comparison with boundary conformal perturbation theory 
(BCPT), boundary truncated conformal space approach (BTCSA) and the boundary 
analogue of the Liischer formula. 
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1 Introduction 



The nonlinear integral equation (NLIE) approach PU E] is a powerful tool for studying finite- 
size effects in the sine-Gordon model with both periodic [2] - [H] and Dirichlet [7JIH] boundary 
conditions. A NLIE has recently been proposed P for the ground state of the sine-Gordon 
model on a finite interval with more general integrable boundary conditions [TlUlllj. including 
new boundary terms proportional to the first time derivative of the field {d y <p). We propose 
here a NLIE for some bulk excited states of this model, using which we numerically compute 
the dimensions of these states as a function of scale (the product of the length of the interval 
and the soliton mass) from ultraviolet (UV) to infrared (IR). We perform checks of the UV 
and IR limits analytically. Other approaches to studying this model (although without the 
dyip boundary terms) have been considered in [T2*]-|15|. 

Our NLIE is based on the Bethe Ansatz solution ^01 113 EE] of the XXZ model with 
general (both diagonal and nondiagonal) boundary terms jTHj. A significant limitation of 
this solution is that the boundary parameters are not all independent, as they must satisfy 
a linear constraint relating the left and right boundary parameters. (Such a constraint does 
not arise in the case of diagonal boundary terms [201 EH E2j-) Consequently, our NLIE 
is applicable only when the boundary parameters of the sine-Gordon model (including the 
coefficients of the d y <p boundary terms) obey a corresponding constraint. 

Three different sets of boundary parameters are introduced in the course of this paper: 
the UV parameters (/x± , (p Q , k±) appearing in the boundary sine-Gordon action; the IR 
parameters (r)± , i?± ,7±) appearing in the sine-Gordon boundary S matrix; and the lattice 
parameters (a± , b± , c±) appearing in the XXZ spin-chain Hamiltonian. The relations be- 
tween the continuum parameters (fi±,<p^) and (r]±,$±) are known \F2\ |2~3*]. An important 
challenge in our Bethe-Ansatz-based approach is to have the correct relations between the 
lattice and continuum boundary parameters. Such relations were proposed in Pj. The con- 
sistency of the results presented here for the UV and IR limits of excited states provides 
further support for those relations. 

The outline of this article is as follows. In Section El we collect some results about 
the sine-Gordon model on a finite interval which we use later to compare with the NLIE 
results. In particular, we clarify various aspects of the d y (p boundary terms: the periodicity 
of the coefficients k± (J2.ll)) . and the dependence of the UV conformal dimensions (|2.29j) 
and of the boundary S matrices ()2.35|) on these parameters. In Section El we review the 
construction of the counting function for the corresponding light-cone lattice model [5] , and 
the corresponding expression for the Casimir energy (|4.11|) . Moreover, we derive the lattice 
counting equation (|3.14jl . which is valid also for the homogeneous (A = 0) open XXZ spin 
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chain. In Section 0] we present the continuum NLIE (|4.2|) which follows from the lattice 
counting function. For simplicity, we restrict our attention to source contributions from 
holes and special roots. We also note the relations ([4.12)1 . (j4.13|) between the lattice and 
continuum boundary parameters, and the constraints (|4.14jh (j4.15|) that these parameters 
must obey. In Section 14.11 we analyze the UV limit. We give the NLIE result for the UV 
conformal dimensions of states with arbitrary numbers of holes and special roots (|4.2(J|) . 
and show that it can be consistent with the CFT result (|2.29|) for appropriate values of 
the boundary parameters. In Section T4.2I we analyze the IR limit. In particular, we verify 
that the IR limit of the NLIE for a one-hole state is equivalent to the Yang equation for 
a particle on an interval. A noteworthy feature of this computation is that the boundary 
5* matrices JT] which enter the Yang equation are not diagonal. Our numerical results, 
including comparisons with boundary conformal perturbation theory (BCPT), boundary 
truncated conformal space approach (BTCSA) j2U and the boundary Liischer formula 
[26| 127] , are presented in Section 03 Section |H1 contains a brief summary and a list of some 
remaining problems. In Appendix A we present a discussion of BCPT and BTCSA. 

2 The sine-Gordon model on a finite interval 

In this Section, we collect some results about the sine-Gordon model on a finite interval 
which will be needed later for making comparisons with NLIE results. In particular, we 
clarify various aspects of the d y (p boundary terms: the periodicity of the coefficients k±, the 
dependence of the UV conformal dimensions on these parameters, and the dependence of 
the boundary S matrices on these parameters. 

2.1 Action 

Following we consider the sine-Gordon quantum field theory on the finite "spatial" 
interval 1, with Euclidean action 

dy / dx A((p , d^ip) + dy 

-OO JX— J — oo 

where the bulk terms are given by 

A(ip , df.if) = ^{d^) 2 - n bu i k cos{f3ip) , (2.2) 



B_((p,d y <p) _ +B + (<f,dy<f) 



(2.1) 
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and the boundary terms are given by 1 



B±(<p , d y (p) = -/j± cos(^(v3 - <Pq)) ± iK±d y Lp . (2.3) 

As noted in [9], this action is similar to the one considered by Ghoshal and Zamolodchikov 
[TT] . except that now there are two boundaries instead of one, and the boundary action 
()2.3j) contains an additional term depending on the "time" derivative of the field. In the 
one-boundary case, such a term can be eliminated by adding to the bulk action ()2.2j) a term 
proportional to d x d y (p, which has no effect on the classical equations of motion. However, 
in the two-boundary case, one can eliminate in this way only one of the two k± parameters 
(say, k+), which results in a shift of the other ( 

The k± parameters are real. The factor of i in the d y <p terms in ()2.3j) (which was missed 
in [Hj) is introduced by the Wick rotation from Minkowski to Euclidean space. Indeed, the 
Minkowski-space action is given by Am = S^oo^t L M , with 



£ + dx l - ([d tV f - {d xV f) - K+d t cp(x + ) + k.OMx.) + . . . , (2.4) 



where the ellipsis (• • •) represents the mass terms (proportional to fibuik or fi±) which we 
have suppressed for brevity. With ip and k± real, the Minkowski-space action is real, as is 
necessary. Rotating to Euclidean space t = —iy, d t ip = id y <p, we see that 



"+ 1 1 

dx y ((dyV) 2 + (dx<p) 2 ) + iK + d y ip(x + ) - iK-d y (p(x-) + ...> = -L E ■ (2.5) 



Hence, Am = I^oo(~ ^v) (~Le) = iAe, with the Euclidean action Ae = S^oa^V Le given 
by Eqs. (|2~T |l -(j23 |l . As usual, e iAM = c~ Ae . 

An important observation is that the k± parameters are periodic, with periodicity j3/2. 
Indeed, first observe that the action (|2.1|) - (j2.3|) has the periodicity 2 

47T 

<p{x,y)^p(x,y) + — . (2.6) 

The contribution from the d y cp boundary terms to e~ AE in the Euclidean path integral 
/ T>ip e~ AE is evidently given by 



1 While in |H] the coefficients of d y ip are expressed in terms of the parameters y± in the boundary S 
matrices l|2.32|l . here we instead denote these coefficients by new parameters k±. 

2 Although the bulk terms (|2.2|l have the periodicity ip(x,y) i— > (p(x,y) + =§■, the boundary terms l|2.3|) 
have only the reduced periodicity l|2.(j|) . 
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where 

/oo 
dy d y cp(x, y) = tp(x, y = oo) - <p(x, y = -oo) . (2.8) 
-oo 

Let us compactify the y axis to a circle, so that y = — oo and y = oo correspond to the same 
point. It follows that the sine-Gordon field on the boundary at y = — oo must be identified 
with that at y = oo, up to the periodicity ()2.6|) . Hence, 

An 

Ay?(x ± ) = — n± , (2.9) 

where n± are integers. It follows that the contribution (J2.7|) to e~ AE becomes 

e - 4 r (K+n+ - K - n - ) , (2.10) 

which has the periodicity 

k±^ k± + ^. (2.11) 

We recall here that it is useful to introduce the parameters A and v which are related to 
the bulk coupling constant [3, 

8vr 1 , 

A "- 1 = — • (212 > 

Hence, the attractive (0 < (3 2 < 4tc) and repulsive (47r < [3 2 < 8n) regimes correspond to 
the ranges 1 < v < 2 and v > 2, respectively. 



2.2 Ultraviolet limit 

The sine-Gordon model (|2.1|) - (j2.3|) can be regarded as a perturbed c = 1 boundary conformal 
field theory (CFT). In the ultraviolet limit fibuik, H± — > 0, the Minkowski-space Lagrangian 
is given by (see ()2.4|0 

f x + L / \ 
L M = J dx - ((d t ip) 2 - (d x f) 2 j - K+d t <f(x + ) + K_d t ip(xJ) . (2.13) 

It follows from the variational principle that tp(x,t) obeys the massless free field equation 
and Neumann boundary conditions, 

(d 2 - df) <p(x, = 0, d x tp{x, t)\ x=x± = . (2.14) 

Although the dt<p boundary terms do not affect the central charge (they are "marginal" per- 
turbations), they modify the expression for the conformal dimension, which we now proceed 
to compute by canonical quantization. 
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The canonical momentum II conjugate to (p is given by 

6H 

n = 777; — r = 9 t (f - k + 5(x - x+) + K_5(x - xJ) , (2.15) 
6(dt(p) 

where £ is the Lagrange density whose spatial integral is the Lagrangian (|2.13|) . We expand 
if in terms of modes, 

<p{x, t) = Vo + Troy + 4= E -<*n cos(7rn(x - x.)/L) e - i7rnt/L , (2.16) 

where L = x + — X-. One can verify that this expression satisfies the equations of motion 
(I2.14|) . The mode expansion for IT is obtained by substituting (I2.1fij) into (|2.15|1 . Note that 
the momentum zero mode il is given by 

n = f + dx il(x, t) = 7T - K + + K_ . (2.17) 

The canonical equal-time commutation relations 

[U(x, t) , <f(x', t)} = -i5{x - x') , [U(x, t) , U(x', t)} = [tp(x, t) , ip(x', t)]=0, (2.18) 
imply that [a n , a m ] = n5 n+mfl and 

\no,<p ] = -i. (2.19) 

The Hamiltonian is given by 

H= [ X+ dxH, H = Udt(p-C. (2.20) 



Substituting the mode expansions, we obtain 

H = — — 7Tq + modes , (2.21) 

where "modes" represents the contribution of the oscillators a n . The wave functional of the 
zero mode is a plane wave, 

(y> ) = e in °^° . (2.22) 

Let us now compactify the Boson on a circle with radius r, which means that the theory 
is invariant under 

<p(x, t) i-> <p(x, t) + 2irr , (2.23) 
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or equivalently, ipo i— > ipo + 27rr. Imposing this condition on the wave functional ^(tpo 
^(<Po + 27cr) implies the quantization of the momentum zero mode, 

n 



n 



(2.24) 



where n is an integer. In view of (j2.17|) . (|2.21j) and ([2.24)1 . the zero-mode contribution to the 
energy is 



1 (n 

E n = — \ - + - K,- 



2L \r 

Comparing this result with the CFT result 

7T , 



7T 



E 0,n = -^jf&ff ~ l ) = £ Afl 

leads to the following expression for the conformal dimension 

1 fn ^ 2 

&n = — [- + K+ - K. 



(2.25) 



(2.26) 



(2.27) 



For the boundary sine-Gordon model and its UV limit, the compactification radius must 



be 



corresponding to the periodicity (|2.6|) . We conclude that A n is given by 



(2.28) 



1 fn/3 
= 2~n [T + K - 



1 



Aviv - 1) 



2n{u - 1) + -{y - 1)(«+ - k_) 

P 



(2.29) 



This result is consistent with the k± periodicity (|2.11j) . Also, this result is "dual" to the 
corresponding result for a free massless Boson with Dirichlet boundary conditions 



2.3 Boundary S matrices 

Results from the theory on the left half line ^T] imply that the right and left boundary S 
matrices R{6 ; T}±, i?±, j±) are given by 

R{0 ■ V , 0, 7 ) = r„(0) ri (9 ; V , 0) M(6 ; V , 0, 7 ) , (2.30) 

where M has matrix elements 

M(9; V ^)=( mn mi2 ), (2.31) 
6 



where 



mn = cos 1] cosh •§ cosh(A6 l ) + i sin 77 sinh $ sinh(A#) , 

rri22 = cos 77 cosh $ cosh(A#) — i sin 77 sinh $ sinh( X6) , 

m.12 = ze n sinh(A#) cosh(A6 l ) , 

m 21 = i e ~^ sinh(Afl) cosh(Afl) . (2.32) 

Moreover, the scalar factors have the integral representations (see, e.g., [H]) 

du sirm((i/-2)o;/2)sinh(3a;/2)l 



roc dfjj 

r (6) = exp<2i / — sin(2#a;/7r) 
Jo U) 



sinh((v - l)w/2)sinh(2w) 



ri(e;v,#) = —a(r),e)a(i#,e), (2.33) 

cos rj cosh v 



where 



f f°° du ,, nS ., .„ .. .. cosh((i/ — l)a;x/7r) 1. 

0) = exp \ 2 / — sin ((in - 6 W 2tt sin (0a;/ 2vr — — — v \ ; ( — / , , 1(2.34) 
v ' ; F | io w u ; /v v /v ;; s i n h((z/ - 1)cj/2) cosh(u;/2) J v ; 

Note the presence of the factors e ±J7 in the off-diagonal matrix elements m 12 and 77121, 
which are related to the presence of the d y tp terms in the boundary action ()2.3|) . and which 
are absent in the case of a single boundary jTT]. In [5| an argument from JT] was borrowed 
to determine the relation between the (real) parameters j± in the boundary S matrix and 
the (real) parameters k± in the boundary action; namely (after correcting for the missing 
i), 7± = /3k±/tt. However, this relation seems to be incorrect, since it would imply that the 
factors e n± do not have the periodicity ([2. II)) . We propose here instead the relation 3 

An 

7± = -J K ± > (2-35) 

which implies that the factors e* 7± (and hence, the boundary S matrix) do have the expected 
periodicity ()2.11|) . 

The relation ()2.35|) for the right boundary can be understood from elementary consider- 
ations. Indeed, for the theory on the left half-line x_ = —00 , x + — 0, the Minkowski-space 
Lagrangian (j2.13|) can be written in the form 

L M = L M (k+ = 0) - K+d t (p{x = 0) . (2.36) 

The amplitude for a process can be expressed using a path integral of the form 

J Vtp exp (i J°° dt L M \ = e - iK+ALp J Vcp exp (i J°° dt L M (k + = 0)) , (2.37) 

3 The two relations coincide at the free Fermion point. 
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where now 

Aip = <f(x = 0, t = oo) — (p{x — 0, t — — oo) , (2.38) 

with appropriate initial and final configurations of the field (p at t — =Foo. For definiteness, 
we can fix the asymptotic condition (p(x = — oo) = 0. For a process involving a soliton 
reflecting back into a soliton, we have 

(p(x = 0,t = Too) = — => A^ = 0. (2.39) 

P 

Hence, for such processes the amplitude is independent of k + ; and this is also true for the 
reflection of an antisoliton into an antisoliton. For a soliton reflecting into an antisoliton we 
have 

27T An 
V (x = 0,t = T°o) = ±j => A V = -j, (2.40) 

which results in a phase factor 

exp • (2-41) 

For an antisoliton reflecting into a soliton the resulting phase factor is the inverse of the 
above. This leads to the relation ()2.35|) between the parameters in the Lagrangian and the 
reflection factor for the right boundary. 

Similarly, for the left boundary, we consider the theory on the right half-line x_ = 
0,a; + = oo. The corresponding Lagrangian and path integral are given by ()2.36j) and ()2.37jl 
with k + i — > —K-] and we now fix cp(x = oo) = 0. For a soliton reflecting into an antisoliton 
we have 

27T 47T 

^( x = 0,t = T oo) = Ty A( ^ = ^"' ( 2 - 42 ) 

which results in a phase factor 

exp , (2-43) 

and leads to the relation (j2.35j) for the left boundary. 

We find further support for relation (|2.35|) from a study of the UV and IR limits of the 
NLIE in Sections IP and 

The relation of the boundary S- matrix parameters r)±,'&± to the parameters n± , <Pq in 
the boundary action (|2.3j) is given by O 

aoa[—{T}± + v&±)\ = ^e^° , 

cos | — (n± - itf±) ] = ^e^o (2.44) 
\8tt J fi e 
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where 



2{jt 



bid k 



\ * (f ) ' 



(2.45) 



Note that we have introduced an additional minus sign on one of the boundaries. That is, 
the UV-IR relation is different on the two boundaries, the difference being in the sign of 
<p . The two sets of boundary parameters (fi± , y?o , and (t)± , i?± , 7±) can be regarded 
as "UV" and "IR" boundary parameters, respectively; hence, the relations ()2.35|) . (|2.44|) . 
()2.45|) correspond to UV-IR relations. 



3 The lattice counting function 

The light-cone lattice |3TJ1 l3*H l3*2*j version of the sine-Gordon model is similar to the XXZ spin 
chain, the main difference being the introduction of an alternating inhomogeneity parameter 
±A. The solution [TBI E! leads to the Bethe Ansatz equations [Oj 

/»M(A,-) = if\ j = l,...,M( + ), (3.1) 

where } are integers, and the lattice counting function h^ + \X) is given by 4 

h (+ \\) = ^{N[q 1 {X + A) + q 1 (\-A)] + q 1 {X)+r 1 {X) + q 2a __ 1 (\)-r 1+2ib _(X) 

+ g 2a+ -i(A) - r 1+2lb+ (\) - fe(A - A fe ) + q 2 (X + A fe )] } . (3.2) 



The functions q n (X) and r n (A) are odd, and are defined by 



k=l 

5 



sinh(f(A+f) | 

q n (X) = 7r + zln Y~ (- = 2 tan (cot(n7r/(2z/)) tanh(7rA/V)) 

sinh(£(A-f; 



cosh(^(A + f) 

r„(A) = 2 In ^ (3.3) 

ccsh(f(A-f)) 

The real lattice boundary parameters a±,b±,c± must satisfy the constraints 

a_ + a + = ±|c_ — c + | + , 

+ = 0, (3.4) 

4 It should be clear from the context whether the symbol A refers to the value (|2.12|) of the bulk coupling 
constant or to the rapidity variable, as in (|3.2J) . 

5 The branch cut of lnz is chosen along the positive real axis; hence, ln(— 1) = in. 
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where the integer k G [—{N + 1) , N + 1] is even if N is odd, and is odd if N is even. The 
parameters a± can be restricted to the fundamental domain \2a± — 1| < 1v. The number 
M^ + ' of Bethe roots is given by 

M^ = ^(N-l + k), (3.5) 

where k is the integer in the constraint (|3.4j) . 
The corresponding energy is given by [7J E2] 6 

£ = ~t E [ai(Ai + A) + oiCAj - A)] , (3.6) 

° 3=1 

where 5 is the lattice spacing, and 

Id ... 1 sin(n7r/V) 

«n A = ^TT^(A) = \ \ . (3.7) 

27r aA z/ cosh(27rA jv) — cos(mr / z/j 

For given values of the bulk and boundary parameters, the counting function h^ + \X) does 
not give all 2 N energy levels. The remaining levels can be obtained from a corresponding 
counting function h^~'(X) with the boundary parameters negated, 

(a± , 6±) (-a± , -6±) , (3.8) 

and with the number of Bethe roots equal to = | (N — 1 — k) PH] • We shall refer to this 
other counting function as the "negated" counting function. Although below we generally 
explicitly discuss only h^^X), corresponding results hold also for h('(\). 

Since the counting function ()3.2|) is odd and has the periodicity A i— ► A + iv, we can 
restrict the Bethe roots Xj to the following region of the complex A plane jH] 

j^eAX), ~ < 3m A < ^| [J |5Re A = , < 3m A < ^| . (3.9) 

The origin (A = 0) is excluded since the corresponding Bethe state would vanish. 

The summation over all the roots in the counting function involves the function 92 (A), 
whose fundamental analyticity strip is |3m A| < min(l ,u — l). Hence, it is useful to classify 
the Bethe roots Xj in the region ()3.9|) as either real, "close" (0 < |9fm Xj\ < min(l , v — 1)), 

6 Thcre is a misprint in the formula (2.29) in |5] for the boundary energy of the (homogeneous) open XXZ 
chain: the first term in the second line should be 

^ sinh((^-|2 Q± -l|V/2) 

sgn 2a± - 1 — — , 

smh(jAj/2) 

where sgn(n) is the function defined in i|3.f f \ . 
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or "wide" (min(l , v — 1) < \$sm Xj\ < |). Real solutions of h^ + '(X) = integer which are 
not Bethe roots are called "holes". If an "object" (either a root or a hole) has rapidity 
Xj for which the counting function is decreasing (-^h^ + \Xj) < 0), then the object is called 
"special". We denote by M R , M c , M w , N H and N s the number of real roots, close roots, 
wide roots, holes, and special objects, respectively. Note also that h^^X) is continuous on 
the real A axis. For further discussion about general properties of the counting function and 
the classification of roots and holes, see e.g. [HE]- 

We now proceed to derive a so-called lattice counting equation, which relates M c , M w , 
Nh and N$ (but which is independent of M R and N) for any Bethe state. To this end, we 
first compute the asymptotic limit of the counting function, and take its integer part, 

L/i (+) (oo)J = M (+) + 1 + -(s+ + s^)-k + sgn(z/ - 2)M W + - i (o+ + a_ - k)\ ,(3.10) 

where |_ J denotes integer part, and s± = sgn(a± — |), where the sign function sgn(n) is 
defined as 

!— ■ u^O 
n n ■ (3 - U) 

: n = 

In obtaining the result (|3.10p . we have used the facts 

mi 

q n {oo) = sgn(n)7r for < \n\ < 2u , 



v 



r n (oo) = , (3.12) 

v 

as well as the relation ()3.5)1 to eliminate N, and the second constraint in (|3.4)l . On the other 
hand, one can argue that (see e.g. jl]) 

L/i (+) (oo)J = Nfj + M R - 2N S . (3.13) 

Using the evident relation M (+) = M R + M C + M w to eliminate M R on the RHS of (JHZE3J), 
and then combining with (|3.10j) . we finally obtain the lattice counting equation, 

N H - 2N S = M C + 2M w step(u - 2) + 1 + -(s+ + s_) - k + - - (o+ + a_ - fc)J (3.14) 

where the step function step(n) is defined as 

\ 1 = n > , 
step(n) = " (3.15) 

10 : n < 

The lattice counting equation ()3.14|) is valid also for the homogeneous (A = 0) open XXZ 
spin chain. 
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As a simple example of the utility of this result, consider (as in |U|) the case that N is 
even with k — 1, and look for purely real solutions with no holes or special roots (Mc = 
My/ = Nh = Ns = 0). The lattice counting equation implies 

= + s -) + l\~l K + a- " 1)J > (3-16) 

which is a condition on the boundary parameters a± that is necessary for such solutions to 
exist. Numerical checks suggest that this might also be a sufficient condition for the existence 
of such solutions. 



4 Nonlinear integral equation 

The lattice NLIE can be derived from the lattice counting function ()3.2|) by standard ma- 
nipulations - jHl- The continuum limit consists of taking the number of spins N — > oo, 
the lattice spacing 6 — > 0, and the inhomogeneity parameter A — > oo, in such a way that the 
length L = x + — x_ and the soliton mass m (whose relation to fibuik is given by (|A.12|1 ) are 
given by 

L = N5, m = - e - 7rA , (4.1) 
o 

respectively. Changing to the rescaled rapidity variable 8 = 7rA, and setting p + \8) = 
2mh^ + \8), one arrives at the continuum NLIE for f^ + '(8) 

= 2tmL sinh 8 + iP^ y (8) +ig(0) 

+ til d& Sm G{0-6'-ie) ln(l - e /(+)(e ' +ie) ) , (4.2) 

7T J-oo 

where G(8) is given by 

1 r°° 

G(6) = — duo e- iuj9/7T G(u) , (4.3) 



2lX J-oo 

and the Fourier transform G{uj) is given by 

C(u) = sinh((z, - 2V/2) 

y) 2sinh((z/-l)^/2)cosh(cu/2) ' K ' J 

Furthermore, P£Jy(9) is the odd function satisfying P^{9) = 2R^(9), where i? (+) (#) is 
given (as in ([4.3)1 ) in terms of its Fourier transform 

£(+) ^ = sinh((z/- 2)^/4) cosh(^/4) smln> - 2)u/2) 

{ ' sinh((v- l)w/2)cosh(w/2) 2 sinh((z/ - l)w/2) cosh(o;/2) K '' 

s + sinh((i/- [2a + - l|)o;/2) sinh((l + 2i& + )u/2) 

+ 2sinh((i/- l)u;/2)cosh(u;/2) + 2 sinh((z/ - l)u;/2) cosh(u;/2) + ^+ ^ ~) ' 
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where (+ <-» — ) is a shorthand for two additional terms which are the same as those on the 
second line of (|4.5|) . but with a+, s + and 6+ replaced by a_, s_ and 6_, respectively. 

Moreover, <7(0) is the source term. For simplicity, we henceforth restrict our attention to 
source contributions from holes and special roots; other bulk sources (close or wide roots) 
can presumably be treated in the same manner as in |H El- The source term is therefore 
given by 



9(0) = £ - *f ) + X(* + *f )] (4-6) 

3=1 
N S 



where x(0) is the odd function satisfying x'(0) = 2G(6 I ), the latter function being given by 
(|4.3j) . Finally, Of and Oj are the positions of the holes and special roots, respectively, whose 



corresponding distinct, positive integers we label by if and if, 

fM{9f) = 2mlf , fW{0f) = 2mlf . (4.7) 

For the continuum model, the value Ij = is excluded because the corresponding rapidity 
9j = is not physical. 

The energy is given by 

E tbulkL -\- ^boundary ^Casimir i (4-8) 

where e buik and e boundary are given by 

tbuik = ^m 2 cot(z/7r/2) (4.9) 



and 



m 

^boundary ~ 



2 

and Ecasimir (order 1/L) is given by 



/ i , «. , cos((f — 2s_i_a_i_)7r/2) cosh(7r6 + ) . , 
-cot (vir 4) - 1 + u , +) 1 > + . / + + (+<-- 

sm{u7r/2) sm(z/7r/zj 



,(4.10) 



11 H iVg 

Ecasimir = TU £ COsh Of - 777, £ [cosh^ 5 + it) + COsh(6>j 
3=1 3=1 

— f°° rift C^rr, cinVi ( ft J- icl lnH _ c / C+) (0+^ 



— / d0 Qfm sinh(fl + it) ln(l - e /v ^ +teJ ) . (4.11) 
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The following relations between the continuum IR (r)± , , d±) and lattice (a± , b±) boundary 
parameters were obtained in [Hj 7 

ti{s±v — 2a±) 



V± = =F- 



2(1,-1) 



*± = (4.12) 

These relations were obtained by comparing the continuum f2J E3] and lattice expressions 
for the boundary energy. A relation between / y± and c± was also conjectured in jU] 

7 ± = ^, (4-13) 

for which we shall find further support from analysis of the UV and IR limits in Sections 14.11 
and 14 .21 below. Note that these relations together with the constraints (J3.4j) among the lattice 
boundary parameters imply corresponding constraints among the continuum IR boundary 
parameters, 



77- -77+ = T|7--7+l + 



^( s + + s-)v - k 



v - 1 

= 0. (4.14) 

The UV and IR limits also imply that, in the continuum limit, the integer k must be restricted 
to odd values, 

k = odd . (4.15) 

The UV-IR relations ()2.35|) . ()2.44j) . ()2.45|) imply corresponding constraints among the UV 
parameters appearing in the boundary sine-Gordon action. We emphasize that our NLIE 
describes the sine-Gordon model only for values of boundary parameters which satisfy these 
constraints. 

It is interesting to consider the continuum version of the lattice counting equation ()3.14j) . 
For the case of periodic boundary conditions, Destri and de Vega have argued [I] that 
the continuum result coincides with the lattice result without the integer part term. It is 
therefore plausible that for the case at hand the continuum counting equation is given by 

N H - 2N S = M C + 2M w step(u - 2) + 1 + ^(s+ + s_) - k . (4.16) 

We shall find some support for this conjecture when we consider the UV limit below. 8 

7 We have already implicitly noted the relation between the continuum and lattice bulk parameters in 
(|2.12|l . Indeed, in 0, v is defined as v = n/fi, where /j£ (0 , 71") is the bulk anisotropy parameter of the XXZ 
spin chain. Here we do not explicitly introduce the lattice anisotropy parameter /i in an effort to reduce the 
number of parameters appearing in the paper. 

8 We recall 0] that special objects cannot appear in the IR limit. They can appear as mL — ► 0, in such a 
way that Nn,eff — Nh — 2Ns remains constant. 
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4.1 Ultraviolet limit 



We now consider the UV limit mL — > 0. In this limit, only large values of \8\ contribute to 
the Casimir energy. 9 Indeed, the driving term of the NLIE implies that one must consider 
the scaling limit 

d = ln^-+6, (4.17) 
mL 

where is finite as mL — > 0. Hence, one must distinguish whether the positions of the 
sources remain finite or scale in the same manner. Let (Nf^) be the number of holes 
whose positions Of {Of °°) remain finite (scale), with corresponding integers If {If 00 , 
respectively); and similarly for the special roots. Hence, 

of ™ = \n^- + 9? 00 , ef 00 = ln^- + ef 00 , (4.18) 

J mL J J mL J 

with Of 00 and Of 00 finite as mL — > 0, and 

N H = N° H + N™, N S = N° S + N™. (4.19) 

Proceeding to compute the Casimir energy from (J4.2)) and (J4.ll)) as in @]-[7], imposing the 
constraints (|3.4jh and recalling that Ecasimir = w ^h c e // = 1 — 24 A n , we obtain 

-{s+ + s_) + iV£ - - 2N% + 2iY|° - (jfe - 1 + 2N° H - 4iV° 



!2 A 1 "^ 
+ E J f °° - 2 E J / °° - - 2N^)(N% - 2N™ + 1) . (4.20) 

i=i i=i 2 

The result (J4.20)) for the conformal dimension is consistent with the CFT result ()2.29j) if 
we identify 

c±=~(u-1)k±, (4.21) 

and set 

^(s+ + s_) + JV& - iVf - 2iV£ + 2iVJ° = fc - 1 + 2iV° - 4iVg = 2n , (4.22) 

where n is an integer. The relation (|4.21|) together with ()2.35|) implies the relation ()4.13|) 
between the boundary parameters j± and c±. 10 Note that the result ()4.20j) for the conformal 

9 Since the counting function is odd, we consider explicitly only 9 — > 00, and we double the result in order 
to also account for the # — > — 00 contribution. 

10 We remark that H4.21| and the periodicity 1)2.11(1 imply the periodicity c± — ► c± + 2(v — 1); i.e., the 
periodicity e± — > c± + 2v of the lattice model |0] becomes "renormalized" in the continuum. 
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dimension contains not only the contribution ()2.29|) corresponding to the highest weights, 
but also additional integer-valued terms corresponding to their descendants. 

It follows from (ETTHD and (ET22J) that 

N H -2N s = l + ^(s + + s.)-k, (4.23) 

which is in agreement with the conjectured continuum counting equation (J4.16)) for the 
special case Mc = Mw = that we are considering. Moreover, ()4.22|) implies the result 
(j4.15|) that k must be an odd integer. 

Here are three examples: 

Ex. 1 A Bethe state consisting of only real roots and no holes or special roots (Nh = 
Ns = 0) with k = 1 has, according to ()4.23|) . a "good" UV limit only if the boundary 
parameters satisfy s + + s_ = 0; and its UV dimension, according to ()4.22|) . is given by 
()2~2T$j) with n = 0. 

Ex. 2 A Bethe state consisting of only real roots and 1 hole and no special roots (Nh = 
1 , Ns = 0) with k = 1 has a "good" UV limit only if s + + s_ = 2. Its UV dimension 
is given by ()2.29j) with either n = 1 or n = 0, depending on whether in the UV limit 
the hole's rapidity remains finite (N H = 1 , N^ = 0) or infinite (N H = , Nf? = 1), 
respectively. 

Ex. 3 For k — — 1, it is possible to have a Bethe state with a "good" UV limit consisting 
of only real roots and 2 holes and no special roots (Nh = 2,N$ = 0), provided 
s + + s_ = 0. Its UV dimension is given by (|2.29j) with n values equal to either 1,0, 
or —1 depending on whether in the UV limit the hole rapidities are either both finite, 
one finite and one infinite, or both infinite, respectively. 



4.2 Infrared limit 

We verify in this Section that the IR limit of the NLIE for a one-hole state is equivalent 
to the Yang equation for a particle on an interval. Indeed, in the IR limit mL — > oo, the 
integral terms in the NLIE ()4.2)1 and in the energy formula 1)4.11)1 are of order 0(e~ mL ) and 
can therefore be neglected. For a single hole with rapidity 8h, the NLIE becomes 

f (+ \9) = 2imL sinh 9 + iP^ y {9) + i X {9 - 9 H ) + i X {9 + 9 H ) . (4.24) 

Noting that e^ (+)( - 6,H < 1 = 1 on account of Eq. ()4.7)1 . we obtain the following relation for 9h 



16 



A similar relation can be derived from the negated counting function, 



^Zim L smh tt H bdry^"" I ' Ij _^ (4 26) 

where P^} y {ff) differs from P^) y {9) by the negation of the boundary parameters ()3.8j) . These 
relations should be equivalent to the Yang equation for a particle on an interval of length L, 

e 2imL ^ 9H R(9 H ;r ]+ ^ + , 1+ )R(9 H ;r ] ^^,^)\9 H ,(±)) = \9 H ,(±)), (4.27) 

where the boundary 5 matrices R(9 ; rj±, i?±, t±) are given by ()2.3()j) . and \9h, (±)) denote 
the two possible one-particle states. 11 

In other words, dropping the subscript H of the hole rapidity, the expressions e ( bdry ^ ^ +X ^ ^ 
should be equal to the two eigenvalues of the Yang matrix Y(9), which is defined by 

Y{9) = R{9 ■ ri + , 1+ )R{9 ■ v _, 7 _) . (4.28) 

Indeed, for a state with real roots and one hole, Eq. (|4.16|) implies that k = |(s+ + s_). 
Since k must be odd, it follows that the only two possibilities are k — s± = 1 or k — s± — — 1. 
For definiteness, we consider the former case, k — s± — 1, and therefore, \ < a± < \ 
From the definitions of P^ y (9) and x($)> ^ follows that 



V. 



PbtrM + X(20) = 2 I" — \sm{u9/7c)R {± \uj) + sm{2u9/7t)G{uj)} . (4.29) 
y 7o a; L J 

With the help of the expression ()4.5|) for i?( + )(co>) (and a similar expression with the boundary 
parameters negated ()3.8j) for i?' '(a;)) and the identity 



wah((v - 2)w/4) cosh(i/a;/4) sinh((z/ - 2)w/2) 



sinh((v - l>/2) cosh(a;/2) 2 sinh((i/ - l)w/2) cosh(c;/2) 
2sinh((z/ - 2)w/4) sinh(3u;/4) 



sinh((z/ — l)a>/4) sinha; 



G(u/2), (4.30) 



we obtain 



= exp « f - { 2 3inh '^ ~ ^ (4.31) 

jo a; L sinh((z/ — l)a;/2) sinh(2u;) 

+ sinh((l ± (u - 2a+))u) + sinh((l ± 2ib + )cj) + ^ j 
2 sinh((z/ — l)u;) coshu; 2 sinh((z/ — l)u;) coshu; -* 



11 In the case of Dirichlct boundary conditions, !##,(+)) and |0fj,(— )) would correspond to one-soliton 
and one-antisoliton states, respectively. For this case, a similar approach for computing boundary S matrices 



was considered in 
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Finally, using the relations (j4.12|) between the lattice and continuum IR boundary parame- 
ters, we obtain 

.(^>(„ + x( 2 .)) = ^^ { 2sinh((,- 2)^/2)^/2) 

P io w v 1 sinh ((u - 1 )u/2 sinh (2co) v ; 



+ 



sinh((z/- l)cj/2)sinh(2cu) 
sinh((l =f (f - 1)2t] + /tt)lo) sinh((l ± (z/ - 1)277_/tt)u;) 



2 sinh((y — l)u;) coshu; 2 sinh((z/ — 1)c<j) cosha; 

sinh((l ±{y- l)2id + / sinh((l ± (v - 1)2m?_/tt)u;) i 
2 sinh((z/ — 1)uj) cosha; 2 sinh((z/ — l)a;) coshu; J 

We now turn to the computation of the eigenvalues of the Yang matrix ([4.28)1 . Recalling 
(12.30)1 . we see that the eigenvalues y^ 1 of Y{9) are given by 



V 



(±) 



r {6) 2 n{6 ; V+ , #+)n(0 ; r/_, tf_)A (±) , (4.33) 



where A^ denote the eigenvalues of M{9 ; rj +1 7 + )M(# ; r^_, 7_). Although the ex- 
pressions for are generally very complicated, a remarkable simplification occurs if the 
boundary parameters satisfy the constraints (j4.14j) . (|4.15|) . Indeed, in that case, the eigen- 
values are factorizable into a product of trigonometric functions, 

A (±) 

= cos(-r? + =f i\9) cos(r?_ =p i\6) cos(id + =F i\6) cos(m?_ =F i\9) . (4.34) 
Recalling the expressions ()2.33)) for r (9) and r 1 (6 l ; rj, i?), using the identity 

cos(x =F zA0) a(x, 9) = exp<2i — sm(20o;/7r) , .. v — ' - ' ' \ , (4.35) 

cosx [ Jo to 2 sinh((z/ — 1)UJ) coshcu J 

and comparing with ()4.33jl . we obtain the desired result 

e i P ^y {9)+xm ) = r (9) 2 ri (9 ; V+ , 0+)n(0 ; t?_, tf_)A^ = y (±) . (4.36) 

That is, we have verified that the IR limit of the NLIE for a one-hole state is equivalent to 
the Yang equation for a particle on an interval. We stress that the boundary S matrices 
entering the Yang equation are not diagonal. 

We remark that, for the case of Nh holes, the Casimir energy in the IR limit becomes 

Ecasimir -> mN H . (4.37) 

Indeed, as already noted, the integral term in the energy formula (j4.11)) can be neglected; 
thus, only the first term of that formula survives. Moreover, the hole rapidities go as 9j ~ ^ 
(since mL sinh 9f ~ 1) for L — > oo, which leads to the result ()4.37j) . 
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5 Numerical results 



The NLIE f!4.2j) can be solved numerically by iteration, and the corresponding Casimir energy 
can then be evaluated with ([4. lip . It is not evident how to best present such numerical 
results for the full range of L e (0,oo). The difficulty is that, in the UV limit (L — > 0), 
Ecasimir = ~ ^II diverges and c e // is finite; while in the IR limit (L — > oo), the reverse is 
true: Ecasimir is finite (!4.H7j) and c e // diverges (if the number of holes is not zero). That is, 
neither Ec as imir nor c e// remain finite over the full range of L. Following [32], we consider 
the dimensionless quantity ("normalized energy") 

ir + mL __ 24(7r + mL)' ^ ' ' 

whose UV and IR limits are both finite: 

£ -> A n"^ fOT L ^ ' ( 5 - 2 ) 

£ -> N H for L -> oo . (5.3) 

We have plotted £ as a function of In Z, where / = mL is the dimensionless scale parameter, 
for various states. 12 



5.1 Ground state 

NLIE results for the ground state (0 holes), whose UV limit is discussed in Ex. 1 at the end 
of Section |4~TI are presented in Fig. [T] As expected, the value of — c e ///24 in the IR limit is 
0; and in the UV limit agrees well with the analytical result for Ao given by (|2.29)l . (|4.20jl . 

We define three regions of I = mL in which we further test, with different methods, the 
ground state energy level obtained by numerically solving the NLIE: 

• The UV region is the small volume region, I < 10 _1 ; here we compare it with boundary 
conformal perturbation theory (BCPT). 

• In the intermediate region, I ~ 1, we test it against truncated conformal space approach 

(tcsa) [211123. 

• In the IR region, where the volume is large, / > 10, we compute its Liischer-type 
[2H 123 correction. 

In all regions we obtain a perfect confirmation of the correctness of our NLIE. 

12 For the case of holes, we plot — c e ///24 vs. In/, which also has the limiting values 15.2J1 . (|5.3(l with 
n = N H = 0. 
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5.1.1 UV region 



Combining the formulae from Appendix A which describe the BCPT and NLIE schemes, we 
obtain the small volume expansion of the NLIE ground state energy 

n f 1 

EnLIe(L) = —€bulkL — €boundary + ~f I -Ejo) — ^ 

Note that E\ Q ) is the conformal dimension of the ground state, which is given by ()2.29|) with 
n = 0; that is, E\ ) = f^- Also, as in Sec. IA.31 here A = £^ = The bulk and boundary 
energies are given by Eqs. (jHHj) and ()4.1())1 . respectively. Computing numerically the ground- 
state energy for small volumes, the coefficients c° m 2 ( A_1 ) can be extracted. Tabled shows a 
comparison between the numerically measured coefficients 13 and the exact values calculated 
from BCPT (jA.ll)) . (jA.12|) for various values of the bulk coupling constant v and for the 
same values of boundary parameters used to generate Fig. ^ The agreement is convincing 
and is of the order of our numerical precision. 



V 


NLIE c° m 2 ^- 1 ) 


BCPT c\ m 2 ^- 1 ) 


1.70 


-5.3215286975 


-5.3215288274 


1.80 


-7.4632436186 


-7.4632435914 


1.93 


-19.5148929102 


-19.5148929079 


2.20 


5.6819407377 


5.6819407318 


2.40 


2.4879276564 


2.4879276494 


2.60 


1.4601870563 


1.4601870411 



Table 1: Comparison of NLIE and BCPT results for c° m 2<yA l \ for various values of bulk 
coupling constant v and for boundary parameter values a + = 1.8, a_ = —0.9, and b + = 
-6_ = 0.41444. 



5.1.2 Intermediate region 

In this region, the energy levels are not dominated only by the first few terms in the UV 

expansion; instead, all the higher-order terms contribute the same way. That is, a non- 

perturbative check is necessary. This is provided by a TCSA calculation, which - being 

a variational method - sums up the perturbative series, in which all the coefficients are 

13 Specifically, we computed m^E^LiE for 100 values of I, from I = 10~ 5 to I = 10 -3 , which we fitted to 
the curve l|5.4|l to obtain estimates for c° m^ A ~ 1 ' > . 
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calculated approximately in a finite-dimensional, truncated Hilbert space. The difficulty is 
in the comparison. TCSA works if the dimension of the perturbing operator is small, that is 
when v is close to one. In this domain, however, the NLIE is not convergent. So one has to 
find a proper range, where the NLIE is convergent and the TCSA is reliable enough. In Table 
12] we present results for v = 1.2 and for boundary parameter values a + = 1.2,a_ = —0.2 
and b + = b_ = 0. The dimensionless NLIE ground state energy data are transformed into 
the TCSA scheme by (jA.6|) and are compared to the dimensionless TCSA data for different 
truncation levels and dimensionless volumes, / = mL. 



Volume 


I = 0.7 


I = 0.9 


I = 1.1 


I = 1.3 


m~ l E TC sA with E cut = 10 


-0.32803 


-0.31248 


-0.31030 


-0.31498 


mr l E TC sA with E cut = 12 


-0.32834 


-0.31284 


-0.31072 


-0.31544 


in' 1 E TC s a with E cut = 14 


-0.32857 


-0.31311 


-0.31103 


-0.31579 


m~ l E TC sA with E cut = 16 


-0.32875 


-0.31332 


-0.31127 


-0.31606 


Tn~ (EjSfLIE + CbulkL + €bdry) 


-0.33067 


-0.31559 


-0.31386 


-0.31895 



Table 2: Comparison of NLIE and TCSA results for the ground-state energy, for v = 1.2 
and for boundary parameter values a + = 1.2, a_ = —0.2, and b + = b- = 0. 

We can see that as we increase E cut the TCSA energies approach the NLIE energy from 
above as a consequence of the variational nature of the TCSA. The truncated Hilbert space 
with E cut = 16 contains 6133 states. 



5.1.3 IR region 

Here we check the exponentially small correction to the ground state energy for large but 
finite volumes. This is dominated by the first breather, with mass mi = 2msin(^-), and is 
given by [22] as 

E NLIE (L) = m 1 1+C ° S ^ ~ S ! U f tan 5± tanh ^± tan ^ tanh %e" miL + . . . (5.5) 
V ; 'l-cos^ + sin^ 2A 2A 2A 2A K ' 

In Figure El this correction is checked as a function ofz/=l + l/A and of the boundary 
parameters. On the figure the logarithm of the dimensionless ground state energy is plotted 
against the dimensionless volume. The upper two lines in descending order are the Liischer 
corrections (J5.5j) for v = 1.25, b + = — &_ = 0.1; and with a + = 1.2 , a_ = —0.2 for the first 
line, and a + = 1.1 , a_ = —0.3 for the second line. The lower two lines have parameters 
v = 1.5, a + = 1.2, a_ = —0.2; and with b + = — 6_ = 1 for the first line, while b + = — 6_ =0.1 
for the second line. The various boxes are the data of the numerical solution of the NLIE. 
The agreement is really excellent. 
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5.2 Excited states 

NLIE results for states with 1 and 2 holes, whose UV limits are discussed in Exs. 2 and 3 
at the end of Section 14. 1| are presented in Figs. El and |3J respectively. Note that the values 
of £ in the IR limit are 1 and 2, respectively, in agreement with (|5.3j) . Moreover, the values 
of £ in the UV limit agree well with ()5.2|) and with the analytical results for A n given by 

In particular, for the 1-hole states (Fig. EJ), we consider integer values I H = 1,2,3,4; 
we find that the corresponding hole rapidities 9 H become infinite in the UV limit, and 
thus, all of these states have n = 0. The values of £ in the UV limit are spaced by 1 
on account of the additional integer contribution to A n in ()4.20j) : as I H = I H 00 increases 
by 1, so does A n . Similarly, for the 2-hole states (Fig. 0J), we consider integer values 
(If , If) = (1, 2), (1, 3), (1, 4), (2, 3); we find that both hole rapidities become infinite in the 
UV limit, and thus, the states have n — — 1. As If + If increases by 1, so does the limiting 
UV value of A„. Hence, the lowest line is (lf,lf) = (1,2), the second-lowest line is (1,3), 
and the next two (almost degenerate) are (2,3) and (1,4). 

6 Conclusion 

Starting from the Bethe Ansatz solution P^ITTI ITH] of the XXZ model with general boundary 
terms, we have derived a nonlinear integral equation for some bulk excited states of the 
sine-Gordon model on a finite interval with general integrable boundary interactions ^01 
ITT] , including boundary terms proportional to d y <p. We have used this NLIE to compute 
numerically the dimensions of these states as a function of scale, and have checked the 
UV and IR limits analytically. We have also verified that the ground-state NLIE agrees 
well with boundary conformal perturbation theory (BCPT), boundary truncated conformal 
space approach (BTCSA) and the boundary Liischer formula. An advantage of the latter 
approaches is that they are not restricted to values of the boundary parameters that obey the 
constraints ()3.4|) . ()4.14|) . The consistency of the results provides support for the proposed 
relations between the lattice and continuum boundary parameters. 

The result (|2.29|) for the conformal dimensions of a free massless Boson with Neumann 
boundary conditions and d y (p boundary terms, which is "dual" to the corresponding result 
for a massless Boson with Dirichlet boundary conditions [2EJ EH] , may have applications in 
other contexts, such as string theory. 

There are many issues that remain to be addressed. Among these are the proper treat- 
ment of complex (bulk) and imaginary (boundary) sources in the NLIE. While the former 



22 



problem is in principle understood |U |S], the latter problem is still not well understood 
even in the simpler Dirichlet case [HI EE] • Moreover, it would be interesting to extend the 
comparison of NLIE with BCPT and BTCSA also to excited states. 
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A Boundary Conformal Perturbation Theory 

and Boundary Truncated Conformal Space Approach 

Boundary conformal perturbation theory (BCPT) and boundary truncated conformal space 
approach (BTCSA) [23 EH] can be applied if the theory is a relevant perturbation of a 
boundary conformal field theory: 

L = L BC ft + L pert = L BC ft ~ l^buik / $(x,t)dx - - , 

J X- 

where Lbcft is the Lagrangian of the UV limiting boundary conformal field theory, <$>(x,t) 
is a relevant bulk primary field of weights (h, h) and ty± are relevant boundary fields living 
on the left /right boundaries of the strip with weights A±. For simplicity we will suppose 
that h = h = A_ = A + =: A. and put x_ =0; x + = L. 14 



A.l Hamiltonian approach 

We are interested in the spectrum of the Hamiltonian: 

H(L) = Hbcft(L) + fibuik [ L <f>(x,t)dx + M-*-(t) + • 

Jo 

The volume dependence can be obtained by mapping the system to the upper half plane 
(UHP) via z = e l ^^ x+t \ where t = —iy is the Euclidean time. The transformation rules of 

14 We emphasize that Sees. IA.ll and lA.2l are more general than the main body of the text, as they are valid 
for any perturbed boundary conformal held theory. 
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the primary fields are given by: 



' 2 \ A . \ A 



<$>( X}t )=\—zz\ <t>(z,z) ; * ± {t) = \^-z) , (A.l) 



Changing the integration variable to 9 = jx and taking the Hamiltonian at t — 0, we have: 
H(L) = £ (Lo-^+fiw^y*' 1 J*He»,e^)M 

+ 0-(£)V(l) +/*+(£)%+(-!) , 

where L is the spectrum of the boundary conformal field theory with central charge, c. 
The spectrum of H(L) can be calculated at least in two different ways: using perturbative 
(BCPT) and variational methods (BTCSA). 

In the variational method we use, as input, the eigenvectors, \n), of the unperturbed 
(boundary conformal) Hamiltonian. For practical reasons we consider the eigenvectors hav- 
ing energy less then a given value, E cut , and perform the calculation numerically. (Technically 
this means diagonalizing the truncated Hamiltonian). 

Standard perturbation theory gives rise to the following perturbative series for the energy 
level labeled with its unperturbed UV limiting vector |n), 

7T / c ^ /7r\ fe ( A_1 )\ 

En(L) = I i^E ln) -- + j:^ b ui k ^±,A)[^ I j J , (A.2) 

where E\ n ) denotes the conformal energy on the UHP. For the ground state the first few 
terms have the form 

E Q (L) = - \E\ 0) - — + c?(/i ± , A) [-) + c° 2 (fi bulk , n±, A) [-) +.. 



where 



and 



(A.3) 

c?(/i ± , A) = (0| (/,_*_(!) + |0> (A.4) 



c°(^ fe ,/x ± ,A) = n bulk f <0|$(e* e -*)|O)d0 

Jo 



(0|^(il)|n)(n|^(jl)|0) 

+ 1^ 1^ p — p • ( A - 5 ) 

i,j={±} neH ^10) - ^\n) 

The large volume behavior of the ground state energy is 

E (L) = e bu i k L + e boundary + et oundary + 0(e" mL ) as L —> oo . 
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The ground state energy in the NLIE description, however, is normalized differently as 
Eq LIE (L) = 0(e~ mL ) when L — > oo. The correspondence between the two schemes is 

Eq(L) = Eq LIE (L) + 6bulkL + ^boundary + boundary ■ (A-6) 



A. 2 Lagrangian approach 

The evaluation of the second order term in the Hamiltonian perturbation theory (|A.5|) is 
cumbersome, since we have to sum up the various matrix elements. We can avoid this 
calculation by doing Lagrangian perturbation theory instead. We compactify the strip in 
the time-like direction on a circle of radius R and consider the large R limit of the cylinder 
partition function: 



Z{L, R) = Tr{e~ RH{L) ) = e - REo{L) + ... for R 



oo 



Using the functional integral representation for the partition function with the action, S 

J L(t)dt = SbCFT + Sp er t, 

Z(L,R) = J dl$,ty±]e~ s = J dl$,^ ± }e~ SBCFT e~ Spert 
7 Mg ,g ± ]e- g »<^e-*W _ 

— ^ BCFT ~ _ ^BCFT \C 

^BCFT v 



Zbcft ( X! 



v Spert ) 



\n=0 n - I 

We can obtain the first few perturbative corrections to the ground state energy as 
-~limlog(Z(L,.R)) = E (L) = E» CFT (L) + Y: *<0|*<(0)|0) 

rL | roo 

+fibuik (0\$(x,0)\0)dx- - Yl W / (0|* i (0)* J -(*)|0)d*+... , 

JO 2 iJ={±} J-oo 

where the correlators are the connected BCFT correlators. By transforming the various 
expressions onto the upper half plane we obtain: 

e (l) = E* cFT (L)+Y:vi(l) A mm\o) (a.t) 



+ 



71 N 2A— 1 



L 



/*7T I POO 

Jo 2 i ■ Jo 

where z = u + iv = re ld . Comparing the result with equations (|A.3|A.4|A.5|) we can establish 
the correspondence with the Hamiltonian perturbation theory. Clearly the second order term 
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in (jA.5|) is summed up. One can compare this term directly by inserting the resolution of 
the identity 1 = Y^n \ n )( n \ an d using the conformal transformation property of the fields. 

In any BCFT, using the s/2 invariance of the vacuum, |0), the bulk one point function 
on the UHP can be put to the form 



(0|$(e-e-)|0) = -^ , (A.8) 
while the boundary two point function can be brought to the form 



c 



mwnvuM = |i_ (i i): W |2A > ( A - 9 ) 

where the radial ordering is taken into account. The relevant integrals can be written in 
terms of the beta function B(x, y), both for the bulk and for i — j and for i 7^ j as 



/ duu x ~\l -u) y - 1 = / duu x ~\l + u)~ x - y = B(x,y) 
Jo Jo 



T(x + y) 



The first integral converges only for y > thus < A < | is needed. Collecting all terms, 
the c° coefficient is 

0/ M r(i-A)r(|) 

c 2\^bidki A*±5 A) — HbulkCbulk- 



r(i-A) 

, 2 2 N r(A)r(i-2A) r 2 (A) , A , 

-(^c ++ + // 2 _c_-) pV_ A) -^+ c +-f^AT ' (A ' 10) 



where only the coefficients Q,^, q,- are model dependent. 



A. 3 Boundary sine-Gordon theory 

In the sine-Gordon theory the UV limiting BCFT is described by (J2.13|) . the bulk perturba- 
tion is given by 

$(x,t) = ±(Vp(x,t)+V^(x,t)) ; V (z,z) = n{z,z):e i M*> g >: , 
while the boundary by 

1 / -P^o -^o \ -f> I \ 

*±(t) = -fe-^C/|(t)+e^C/_|(t)j ; V (u) =: e l ^ u ' u) : , 

and A = see [37] for the details. Since the bulk, Vp(z,z), and boundary, U§_(u), vertex 
operators change the eigenvalue of ttq by ft and |, respectively, only even c coefficients 
are nonzero in the expansion (|A.2j) . Moreover, the vacuum expectation value of $ is also 
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zero (|A.8j) . thus the leading perturbative contribution comes from the boundary two point 
function part of (|A.7J) . In the boundary sine-Gordon theory with nonzero k := n + — k_ the 
vacuum is not sl2-invariant and thus (jA.9|) has to be modified. In general for a theory with 
a non s^-invariant vacuum one has to compute the four point functions, instead of the two 
point function, and extract the relevant matrix element from them. Alternatively, in our 
case, one can also use the mode expansion of the field (|2.16j) together with the commutation 
relations (|2.18|) to obtain: 



(0|[/ ± £(Z1)C7 £((j W )|0> 



-_4kA 
+ 3 



This modifies (|A.10|) and gives the leading corrections: 



c° 2 {v±,A,k) = -- 



{ , , ,, r(i-2A) / r(A + ^) r(A-^ 



r(i + ^-A) r(i-^-A) 



P, + _,r(A + ^)r(A-^; 

+fi-ti+ cos-(^ -ip ^ 



2 vru T(2A) 



(A.ll) 



Although the derivation of this formula assumes that < A < 1/2, the final result is analytic 
in A (with possible poles). Therefore it has an analytic continuation for A > 1/2, which, 
since the NLIE is also analytic in A, must coincide with the NLIE result. This is confirmed 
by experience with the NLIE for bulk sine-Gordon and bulk supersymmetric sine-Gordon 
models. Using the UV-IR relation ()2.44j) . ()2.45|) and the mass-gap formula (cf. [SHI) 

^=m-» ^^(^feiy"" , (A.12) 
p tiT(1 - A) \ 2T 



,2(1-A) 

where m is the soliton mass, c^jUi, A, k) can be rewritten in terms of the IR parameters. 
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log(i) 



Figure 1: — c e ///24 vs. log/ for ground state (0 holes), with parameter values v = 1.93, a + 
1.8, a_ = -0.9, and b + = = 0.41444. See Ex. 1 in Sec. HH and Sec. O 
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Figure 2: log(m l E) vs. I for the ground state, with various values of the bulk and boundary 
parameters. See Sec. I5.1.31 
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Figure 3: £ vs. log/ for 1-hole states with integer values I H = 1,2,3,4, and with parameter values 
v = 1.93, a + = 1.8, a_ = 1.9, and b + = -6_ = 0.41444. See Ex. 2 in Sec. O and Sec. I5~2l 




Figure 4: £ vs. log / for 2-hole states with integer values (/-f J^) = (1, 2), (1, 3), (1, 4), (2, 3), and with 
parameter values v = 2.13, a + = 1.8, a_ = —0.9, and b + = — 6_ = 0.50357. See Ex. 3 in 
Sec. El and Sec. O 
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